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ABSTRACT 


A two-dimensional circulation control (4.25 inch diameter by 24 inch span) model 
was designed for installation in the Naval Postgraduate School oscillating flow, low speed 
wind tunnel. An adjustable, tangential blowing slot was included in the design to provide 
a circulation control capability using the Coanda flow effect aft of the spanwise slot. Orifice 
locations were defined for obtaining surface static pressures which could be subsequently 
processed to yield section lift and drag coefficients. 

While the model was being fabricated, pressure system calibrations were made to 
determine the dynamic transfer function from a simulation of the model's static pressure 
orifice when connected to the pressure transducer. Similar transfer functions were 
determined for a static pressure probe. The existing data acquisition system was used to 
process the sampled digital data sequence. 

Clear-tunnel flow calibrations were performed in the (2 foot by 2 foot by 18.6 foot 
long) test section at a mid-length test station using a hot-wire and the calibrated static 
pressure probe. An oscillating velocity was superimposed upon a mean free stream motion 
using an existing rotating mechanism of four synchronized shutters at the end of the test 
section. Analytic estimates of both static pressure and velocity perturbations correlated well 
with experimental results for a frequency range of 14 to 40 Hz with a fixed value of shutter 
blockage. A surprising result confirmed by these test was that, although the amplitude of 
the velocity perturbation wave remained relatively constant at 5 percent of the mean free 
stream velocity, the amplitude of the static pressure wave showed a periodic dependence 


upon the frequency of the harmonic blockage source. 
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I. INTRODUCTION 


A. CIRCULATION CONTROL 

The term "Circulation Control" relates to the development of aerodynamic lift over 
a body immersed in a moving airflow using the technique of tangential blowing of a 
turbulent jet, located on the upper surface of the body. This turbulent jet can be achieved 
by discharging fluid from the interior of the body by pressurizing a plenum with the aid 
of a special blower as illustrated in Figure 1. The location of the tangential injection of fluid 
varies with configuration; i.e., mid-chord region for a bluff body such as a circular cylinder 
and near the trailing edge for an airfoil shape. 

The term Coanda effect refers to the flow entrainment resulting from the injection 
of high energy tangential jet flow into the boundary layer with the effect being an 
alteration of the flow separation point or a shift of the rear stagnation point. The flow 
remains attached aft of the injection slot due to the pressure of the thin layer of turbulent 
fluid until a separation occurs at a point downstream of the slot. The location of the 
Coanda jet separation point depends on factors such as local surface curvature, position 
of the jet slot, the nature of the boundary layer, and the mass flow rate of the jet. 


The tangential jet blowing coefficient, Cy is defined by: 
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where: 


m = jet mass flow rate, 
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Figure 1 Circulation Control of a Bluff Body 


V, = jet velocity, 

Q = dynamic pressure, 

S = reference area. 

A physical interpretation may be placed upon the dimensionless blowing coefficient 
C by recognizing that if the jet were pointing vertically downward and perfect flow 
recovery were achieved, C would translate directly into dimensionless lift coefficient C,. 
lf the jet were pointing aft, C would translate into an equivalent dimensionless thrust 
coefficient. 

The use of circulation control in an airfoil application becomes attractive when it is 
realized that the tangential jet velocity, as measured by C, can produce significant lift 
coefficient changes. Lift amplification, as defined by dC,/dC , can attain values of up to 
50 for well designed airfoil installations. Lift amplification by circulation control on bluff 
bodies such as circular cylinders may reach values on the order of 10 with the reduction 
in amplification from that on an airfoil being attributed to the influence of the large 
Separation region on the aft portion of the bluff body. 

The ability to control lift force by modulating C will be different for different 
applications. In a streamlined airfoil scenario, the concept has been investigated for 
helicopter rotor usage and for fixed-wing aircraft high lift without flaps. C modulation on 
a bluff body has been actively considered for developing a side force upon a rotorcraft tail 
boom where the bluff body cross section is immersed in the downwash flow field from the 
rotor. Thus the time dependent behavior of the circulation control aerodynamics for a 
range of conditions, including different values of blowing and jet slot location with respect 
to the free stream, is of particular interest to designers of the NO TAil Rotor (NOTAR™) 


anti-torque tailboom helicopters, shown in Figure 2. 


As envisioned, a helicopter tailboom equipped with a circulation control form of 
tangential jet flow coming from a horizontal slot would be able to develop a lateral force 
(hence an anti-rotor torque yawing moment), providing the tailboom is immersed in the 
downwash flow field from the rotor. Prototype demonstrations of the concept by 
McDonnell Douglas Aircraft Company have been quite successful. The benefits to be 


derived are the removal of a tail rotor with substantial improvement in reliability. 


B. PREVIOUS EXPERIMENTAL INVESTIGATIONS 
1. Steady Flow Research 

The first experimental investigations with circulation control by tangential 
blowing in steady flow were conducted by Dunham [Ref. 1] and Kind [Ref. 2] in the late 
sixties on circular cylinders and elliptical Circulation Control Airfoils (CCA). The majority 
of experimental circulation control work, for the evaluation of the effect of slot height, 
Reynolds number, trailing edge shape and camber, has been conducted at the David Taylor 
Research Center (DTRC). 

Investigations conducted by Walters [Ref. 3] on a cambered elliptical circulation 
control airfoil, indicated that additional lift augmentation could be obtained with pulsed 
blowing. This produced required lift coefficients at reduced injection mass flow. 

In 1975 a circulation control helicopter rotor was constructed by Kaman 
Aerospace Corporation for concept evaluation on the Navy’s H-2 helicopter. For the 
purpose of evaluation of the rotor performance, sinusoidal pressure waves were used for 
cavity pressurization with amplitude ratios of the order of one, and various combinations 
of leading and trailing edge blowing. High lift-to-drag ratio surface pressure distributions 


were reported. 
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Figure 2 Schematic of a Circulation Control Helicopter Tailboom. 


2. Oscillatory Flow Research 

Unsteady boundary layer behavior was experimentally investigated by Miller 
[Ref. 4]. The instantaneous velocity profiles, in oscillatory laminar boundary layers subject 
to adverse pressure gradients, were measured. Results suggested that the separation point 
moved upstream from the steady-state point and occurred at the farthest upstream point 
at which there was reverse flow at some point in the velocity profile throughout the entire 
cycle of oscillation. 

The performance of a circulation control airfoil operating in a steady free stream 
flow with unsteady blowing was investigated experimentally first by Schmidt [Ref. 5]. An 
experimental technique was established which allowed the measurement of unsteady static 
pressures, with the aid of Scanivalves®, and subsequent digital data processing for the 
express purpose of defining aerodynamic transfer functions of circulation control airfoils. 
Among the results of this investigation were that the airfoil behaved in a linear manner, 
the overall airfoil lift transfer function from harmonic circulation control variations had a 
behavior similar to that of a simple pole in classical control theory, and a frequency 


dependent damping moment was identified. 


C. OUTLINE OF THE INVESTIGATION 

In the present investigation a two-dimensional circulation control circular cylinder 
was designed and constructed for installation in the existing oscillating flow, low speed 
wind tunnel. A calibration technique was performed in order to obtain the dynamic 
transfer function for the unsteady pressure data acquisition system. A further objective was 
to investigate the oscillating flow field produced by a set of rotating shutter vanes aft of 


the low speed wind-tunnel test section. The method of attack consisted of direct 


measurement of clear tunnel test section aerodynamic characteristics in an oscillating flow 
field. Comparison was made with analytical predictions of the oscillating flow field. 

The detailed investigation included: 

1. Designand construction of the circulation control model froma circular cylinder 
for Reynolds numbers below the subcritical point. 

2. Pressure data acquisition system calibrations to determine the dynamic transfer 
function between the pressure orifices on the surface of the circulation control cylinder and 
the signal produced by the pressure transducer. 

3. Analytical prediction of the oscillating flow field parameters and propagation 
characteristics of both static pressure and velocity perturbations due to the reaction of 
harmonic blockage at the end of the test section. 

4. A Static pressure and velocity survey of the wind tunnel test section in steady 
and oscillating flow, without the model installed. 

5. Correlation between analytical and experimental results for static pressure and 


velocity oscillation in the test section, and final results for wind tunnel calibration. 


II]. EXPERIMENTAL EQUIPMENT 


A. WIND TUNNEL 

The investigation was conducted in the Naval Postgraduate School, Department of 
Aeronautics and Astronautics laboratories. The experimental information was obtained 
through the use of the low-speed oscillating flow wind tunnel which is an open-circuit 
indraft facility. 

1. Inlet Section 

The tunnel draws ambient air from the outside environment, through high 

solidity screens, from an eight-foot square inlet into a 16:1 contraction cone. The high 
solidity turbulence screens and the 16:1 contraction ratio result in a very smooth steady 
flow through the test section and help maintain freestream turbulence intensities to less 
than 1 percent for the velocities encountered in the present investigation. 

2. Power Section 

The power section of the tunnel is comprised of two Joy Axivane Fans, mounted 

in series aft of the test section. The fan drive consists of two internal, 75 horsepower, 1750 
rpm, electric motors direct connected to each fan. The fan blades are internally adjustable 
through a pitch range of 25 to 55 degrees, providing a wide operating base. Control of the 
mass flow rate and hence the test section velocity is provided by a set of two externally 
Operated, variable inlet vanes, located immediately upstream of each fan. The multileaf- 
type inlet preswirls the air in the direction of the fan rotation in order to reduce the power 
required to operate the wind tunnel. The range of wind-tunnel velocities is from 10 to 250 


feet per second. 


3. Test Section 

The tunnel test section beginning immediately after the contraction cone, as 
shown in Figure 3, is 223 inches long and has a 24-inch X 24-inch square cross section. The 
test section requires extraordinarily stiff construction dictated by the need to minimize 
deflections due to rapid changes in static pressures when the tunnel operates in the 
oscillating flow mode. For this reason continuous pieces of two-inch thick aluminum slabs 
form the upper and lower test section walls. The front side wall consists of three panels, 
constructed of two-inch thick stress relieved Lucite™, secured by twelve bolts each to the 
top and bottom slabs. The panels are hinged and may be raised by hydraulic actuators, 
providing access to the test section. The rear wall is constructed from two-inch thick 
plywood. 

The lower wall provides ports for the purpose of mounting hot-wires probes, 
pitot-static tubes, and other data acquisition probes, although the plywood construction of 
the rear wall permits mounting of the model and internal or external instrumentation. 

Boundary layer velocity profile measurements at steady flow, have shown that 
velocity variation is less than 1 percent of the mean to within three inches of any wall. 

4. Rotating Shutter Valve 

There are in existence several methods which have been developed for the 
purpose of obtaining an oscillating flow. The method used in this wind tunnel assembly 
was first used successfully by Karlsson [Ref. 6] and later by Miller [Ref. 7]. The method 
uses a rotating shutter valve installed at the exit of the test section consisting of several 
shutter vanes. 

The present shutter vanes consist of four one-inch diameter horizontal steel 


shafts, each of which is slotted to accommodate a flat blade of varying width. The scheme 
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Figure 3 Wind Tunnel Test Section 
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Figure 4 Rotating Shutter Valve Assembly 


11 


























shown in Figure 5, forms an array of four butterfly valves which span the end of the test 
section horizontally. A five horsepower, variable-speed electric motor drives an 
intermediate shaft through a reduction gear assembly. The shaft of the bottom vane is 
driven by the intermediate shaft and then each vane drives its immediate neighbor by 
means of timing belts and pulley arrangements in order to keep all four vanes in phase as 
they rotate. 

The intermediate shaft permits a variety of pulley arrangements and a 
frequency range from two to 240 Hz, while the amplitude of oscillation is controlled by 
blade width with a test section closure from 25 to 100 percent. There is also a provision to 
lock the vanes in the horizontal position to allow the tunnel to produce a steady flow. 

In this experiment, three inch wide blades were used allowing a blockage or 
closure ratio from 8.3 to 50 percent. The test section velocity perturbation varied from three 


to eight percent of the local mean freestream velocity for a frequency range of 10 to 48 Hz. 


B. INSTRUMENTATION 
1. Instrumentation For Pressure Measurement 

A Scanivalve® with its own transducer was used to measure the unsteady static 
pressures from the wind tunnel test section and served as the main element in the design 
of the instantaneous pressure acquisition system for the circulation control model. A 
solenoid controller and an odd-even decoder controlled the opening of the Scanivalve® 
ports to the transducer chamber with the capability to send the Scanivalve® port 
information number to a computer program or to receive port change commands from the 


computer. 
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The analog signal from the Scanivalve® transducer was conditioned by a signal 
conditioning amplifier and was then connected to a recorder board SRA-1200™ (channel 
1). A Tektronix™ type 555 dual-beam oscilloscope was used to monitor the pressure signal 
in real-time. The Scanivalve® unit was mounted on a free-standing mobile dolly in order 
to remain isolated from the wind tunnel, thereby minimizing mechanical coupling and 
signal noise during wind tunnel operation. 

2. Instrumentation For Velocity Measurement 

A linearized, constant temperature hot-wire anemometer described by Miller 
[Ref.8] was employed to measure the velocity oscillations of the clean tunnel test section 
during all experimental runs. The signal from the hot-wire was connected to a second 
recorder board SRA-1200™ (channel 2), and to the same oscilloscope mentioned above 
which was used for pressure monitoring. A Singer Ballantine™ RMS volt meter was used 
to measure the RMS of the velocity analog voltage signal. 

The freestream velocity was also measured from a Pitot-Static probe connected 
to a water manometer with plastic tubing of significant length to damp the pressure 
fluctuations due to velocity oscillations. This system was used to establish a known steady 
flow speed in order to calibrate the hot-wire anemometer and for measuring the mean flow 
velocity during the oscillating flow testing. 

3. Data Acquisition Instrumentation 

The pressure and velocity signals were output continuously from the 
scanivalve® pressure transducer and the hot-wire anemometer, respectively, as analog 
voltages. The signals during recording were stored on the two recorder boards SRA-1200 


(channel 1 and channel 2) with 64K capacity each, installed inside a 80286-based personal 
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computer running at a 12 MHz clock speed. The recording boards were of 12-bit resolution, 
capable of differentiating one part in 4096. 

The data acquisition recorder boards were cued by an electric trigger, which 
was tripped mechanically by a microswitch mounted on the rotating shutter valve system. 
The frequency of the rotating shutter valve was continuously measured and displayed by 
a Dynasciences signal counter, while an Analog Devices™ AD 2501 digital thermometer 
displayed freestream temperature in degrees Fahrenheit, measured by a thermocouple just 
upstream of the rotating shutter vanes. The set-up of the instrumentation used in the data 


acquisition is shown in Figure 5. 
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Figure 5 Instrumentation Set-up For Data Acquisition 
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II. CIRCULATION CONTROL MODEL 


A. GENERAL DESCRIPTION 

The model was a prototype circulation control cylinder, designed and constructed 
in its entirety in the Naval Postgraduate School Department of Aeronautics and 
Astronautics model shop. An aluminum circular cylinder, with 4.25 inches outside and 4.0 
inches inside diameter, was modified to accommodate the injection slot structure. 

Special attention was given during the design process to the injection slot 
configuration. In order to adjust the injection slot width, two separate aluminum pieces 
were bolted inside the cylinder cavity to help five pairs of "push-pull" screws to vary the 
distance between the slot lips. The model end sections were capped by two aluminum, flat, 
circular plates. A one-inch diameter line passed through the rear end plates for external 
air supply to provide pressurization to the cylinder cavity. The cavity may be considered 
aS a reservoir operating at stagnation conditions and provides the supply air for leakage 
through the twenty-inch long tailored slot that creates the Coanda sheet which is initially 
attached to the rounded upper surface of the cylinder. Figure 6 depicts the supply air path, 


the injection slot location, and the main parts of the model. 


B. INSTALLATION IN THE TEST SECTION 

The cylinder was designed to span the 24-inch width of the tunnel test section, 
pivoted at the center of its cross section, for support and injection slot location control. The 
term "Slot Angle of Attack" (Slot AoA) denotes the angle between the cylinder diameter 


which is always perpendicular to the tunnel test section center line and the diameter which 
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Figure 6 Circulation Control Circular Cylinder 
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passes through the injection slot. The slot angle of attack was made adjustable with 
reference to a vernier mechanism, permitting the operator to set any desired angle from - 


40 to +40 degrees relative to the upright direction. 


C. PRESSURE ORIFICE LOCATIONS 

Provisions for twenty-four pressure orifices were made approximately at the model 
midspan station with an even distribution chordwise every 15 degrees. Four additional 
pressure taps were distributed chordwise every 90 degrees, left and right of the main 24 
pressure taps array and seven inches apart. The purpose of these eight pressure taps was 
to verify the two-dimensional nature of the flow field. 

Places for measuring the Coanda jet influence along the model surface, immediately 
after the injection slot, included a static pressure array of seven orifices spaced chordwise 
every three degrees. Figure 7 shows a cross-sectional view depicting the location of the 24 
midspan pressure orifices with respect to the injection slot location, and Figure 8 an upper 
surface view of the cylinder depicting the location of the "Coanda" pressure orifices. The 
decreased angular spacing in the "Coanda" pressure orifices region was for the purpose of 


defining the attached flow region of the Coanda jet. 


D. CONNECTION TO THE PRESSURE MEASUREMENT SYSTEM 

The pressure orifices will be connected to the Scanivalve® ports by means of 
identical lengths of stainless steel tubing 0.065" OD and 0.047" ID. This technique provides 
the same dynamic transfer function for the connection from each of the surface static 
pressure orifices to the simple Scanivalve® pressure transducer. It should be noted that in 
this role, the Scanivalve® acts as a pneumatic multiplexing device since a single transducer 


is selectively connected to the various static pressure orifices using remote control. A 
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Figure 8 Pressure Orifices For Coanda Flow Region 
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calibration procedure was applied to the tubing to determine the transfer function 
necessary to convert the pressure signal received at the Scanivalve® transducer to the 
proper amplitude and phase representative of the pressure distribution over the circulation 


control model. The calibration procedure is described in detail in the following chapter. 
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IV. PRESSURE DATA ACQUISITION SYSTEM CALIBRATION 


The transient response of pressure measuring systems is dependent on two factors; 
the response and sensitivity of the pressure sensor transducer, and the dynamic response 
of the pressure transmitting fluid and connecting tubing. In the experiment the total and 
static pressures associated with the oscillating flow were transmitted to the Scanivalve® 
pressure sensor via a pressure probe and/or stainless steel tubes of specific length. 
Connecting tubing and unavoidable cavities in the unsteady pressure measurement system 
introduced losses and phase lags, causing differences between measured and applied 
pressures. 

Because the pressure transducer inside the Scanivalve® provided a large impedance, 
there was no flow in the tube, and the steady-state average pressure at the end of the tube 
equals the pressure at the probe or the pressure at the pressure orifice location. On the 
other hand for time-varying pressures the tubing acts as a transmission line. 

A harmonic pressure signal input may be considered as being introduced at a 


particular model static pressure orifice in the form of: 


Pt)=P,(1+hcoso) , h<i 


The analog signal sensed by the pressure transducer internal to the Scanivalve® following 
transmission of the information from the model surface through the 0.047" ID stainless steel 


tubing, or from the static pressure probe in the tunnel, may be viewed as being of the form: 


PLonoret) = P,[ 1+ |G(w)|h cos(ot + >] , AK 
where: 

ie = mean pressure, 

h = multiplication factor determining the amplitude of the pressure oscillation 
input, 

|G(@) |= dynamic gain magnitude attributed to pressure tubing and internal 
Scanivalve® configuration, and 

d = phase lag associated with |G(o)|. 


To determine the phase lag and the dynamic gain as functions of frequency, a calibration 
procedure was applied to both the stainless steel tubing and the static pressure probe. The 
static pressure probe was used to measure the Static pressure dynamics in the clear tunnel 


test section. 


A. CALIBRATION INSTRUMENTATION SETUP 

A function generator provided a sinusoidal voltage signal which was serially 
connected, through an amplifier, to a resonator depicted in Figure 9. A pressure signal of 
varying frequency was generated by the 75 Watt acoustic driver of the resonator located 
in the center of the cavity. Two pressure ports at the ends of the cavity provided 
comparative signals with an estimated accuracy of one degree in phase angle. The pressure 
was measured immediately at the first port by a reference transducer and from the other 
port the signal was transmitted to the Scanivalve® pressure transducer via the stainless 
steel tube or the static pressure probe. The two analog voltage signals, after they passed 


through similar signal conditioning amplifiers, were digitized, recorded and stored on the 
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two SRA 1200™ recorder boards. A data acquisition / data reduction program, listed in 
Appendix C , written in Quick Basic”, calculated the phase difference and the amplitude 
ratio between the two Signals using a simple discrete data cross-correlation algorithm. The 
Tektronix™ Type 555 Dual Beam Oscilloscope was used to monitor the same pressure 
Signals in real-time, and a digital voltmeter measured the analog output of the pressure 


transducers to a Sensitivity of +/-0.001 volts DC with a range of +/-10.000 volts. 


B. CALIBRATION RESULTS 
1. Static Transducer Calibration 
Before the dynamic calibration of the stainless steel tubing and the static 
pressure probe was conducted, a static transducer calibration was performed to determine 
the sensitivity of both the reference transducer and the Scanivalve® transducer. A water 
manometer was used to establish static pressure in inches of water and a digital voltmeter 
to measure the transducer’s output in millivolts with an accuracy of +/-0.01 millivolt. Both 
transducers displayed linear static calibration curves, shown in Figure 10, with the 
Scanivalve® transducer having a sensitivity of 0.0346 Volts/psf and the reference 
transducer 0.1063 Volts/psf (three times more sensitive than the Scanivalve®). 
2. Calibration Results For Circulation Control Model Tubing 
Four different tube lengths, 28, 26, 24, and 22 inches, were tested in order to 
forma transfer function data-bank necessary to convert the pressure signals received at the 
scanivalve® back to the proper amplitude and phase representation of the static pressures 
on the surface of the circulation control circular cylinder. For the calibration tests the tubes 
were connected to the Scanivalve® leads with one-inch length plastic tubing, which 


configuration must be preserved in the final model pressure orifices connection to the 
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Scanivalve® leads. Calibration curves showing the phase lag and dynamic gain as a 
function of applied frequency are shown on Figures 11 and 12 respectively. The 28-inch 
length of tubing offers more flexibility for Scanivalve® attachment to the model’s pressure 
measurement system. 
3. Calibration Results For Static Pressure Probe 

The probe used in the clear-tunnel test section static pressure measurements was 
also tested with the same instrumentation setup. A special fitting was made to allow 
positioning of the probe in the resonator cavity. The output lead of the probe was 
connected to the same Scanivalve® used above, by means of a 3 inch plastic tubing, a 
configuration that was preserved during the actual test section static pressure 
measurements. The calibration curve for phase lag displayed a linear behavior, as shown 
in Figure 13, while the calibration curve for the dynamic gain showed a gain increase of 
16 percent of the static value at an applied pressure perturbation frequency of 39 Hz 
followed by a smooth decrease in gain as frequency increased to 100 Hz, as shown in 


Figure 14. 
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Figure 9 Instrumentation Set-up For Dynamic Calibration 
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Figure 10 Static Calibration For Pressure Transducers 
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CONFIGURATION : Stainless Stee! Tubing 0.065" OD , 0.047" ID 
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Figure 11 Tubing Calibration Curve For Phase Lag 
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CONFIGURATION : Stainless Steel Tubing 0.065" OD , 0.047* ID 
and 10° Plastic Tube Connection to the Scanivalye Leads, 
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Figure 12 Tubing Calibration Curve For Dynamic Gain 
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STATIC PRESSURE PROBE DYNAMIC CALIBRATION 
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Figure 13 Pressure Probe Calibration Curve For Phase Lag 
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STATIC PRESSURE PROBE DYNAMIC CALIBRATION 





bn 2 0 0 0 © OD HO HM WM W 
FREQUENCY in Hz 


Figure 14 Pressure Probe Calibration Curve For Dynamic Gain 
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IV. CALCULATION OF STATIC PRESSURE AND VELOCITY OSCILLATIONS IN 
THE TUNNEL TEST SECTION 

Consider the inviscid, unsteady flow through the test section of the tunnel. The air 
is generally moving at the speed u (on the order of 100 ft/sec) with properties P (static 
pressure) p (density) and a (local sonic velocity), while the presence of the rotating shutter 
valve superimposes a periodic variation of velocity on the mean flow. The purpose of the 
analysis that follows was to obtain general expressions for the static pressure and the 
velocity at any instant of time and at any point in the test section, downstream of the 
rotating shutter valve, as a function of the rotation frequency. In order to avoid 
mathematical difficulties, it was necessary to make a number of simplifying assumptions. 
First of all the viscosity and thermal conductivity of the gas were neglected. This means 
that all properties of the gas were related through isentropic relations. Then it was 
assumed that the flow was geometrically one-dimensional, implying that all fluid 
properties were uniform over each cross section of the passage, that no changes in the 
cross section took place in the test section, and that the length of the oscillating pressure 


wave was much greater than the tunnel’s characteristic cross-sectional dimension. 


A. EQUATIONS OF MOTION 


The one-dimensional momentum equation: 


oP 


Gu cu il 
— + Y4— = - 
Oo a pa 


combined with the following expression for the speed of sound (energy equation): 


SZ 


| @ constant entropy 


becomes: 


Op , Kpx) _5 

Ot ax 
Op Op _ du 
—+u— + =0 2 
ie oes re (2) 


Now using small perturbation assumptions, the local density and local velocity have 


the following forms: 


P=py+Ap 


u=U,+Au 


Substitution into equations (1) and (2) of the above expressions for local density and 


velocity, yields: 
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SEP +(UyrAu) OP +(pyrAp)X-0 (4) 





The Taylor series expansion of sonic velocity squared a’, about point p,, gives: 
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Substitution into equation (3) yields: 


a hn ae - 








ay Ap 7 (5) 
Performing an order-of-magnitude analysis because: 
Ap«<p, A  Au«a, 
the following terms: 
auto, A poe no, Pp o-AP~P, , Aum? 20 , Ap <0 


can be neglected and the equations (5) and (4) finally become: 


dAu | OAu 90. OAp _ -0 
al laa Py & 


(6) 


dAp dAp OAu _ 
“a ae PO ae 


Now the perturbations of density and velocity can be written as functions of position 


x and time t: 
A p=p,s(x,t) 


Au=U,r(x,0 


thus the local density and local velocity become: 
p= Poll +5(x,D] 


u=Up[1+r(x,0)] 


Using the above expressions for the perturbed quantities, equations (6) and (7) 


become: 


U Or(x,f) +U 2 Or(x,t) +," ds(x,f) -0 
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and because: 


the above equations take the form: 
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and finally: 


Or(x,t) Or(x,t),  As(x,t) 
mat) +M, — on =0 (8) 


pSsG) , +Mo——— s(x, Gay +My— aoe =0 (9) 
HApf) ax 
Now the following transformation: 
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applied to the equations (8) and (9) yields: 
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Taking the first derivative with respect to x of equation (10), and then the first 


derivative with respect to z of equation (11), yields: 


Q.ar. as 
Re (ee 12 
0a ae? ax? ve 


Os | 


a? 


My<(<)=0 (13) 


The result of the subtraction between equations (12) and (13) is a form of the wave 


equation: 


&s _@s 
ax? a? 


with the well known solution [Ref. 8] : 


S(x,2)=F(x+z) +G(x-z) 
OF: 


= pile + + mele = 
S(z,t) = Fix(1 M af] + G[x(1 M. a,f] (14) 


a7 


Now taking the first derivative with respect to z of equation (10), and then the first 


derivative with respect to x of equation (11), yields: 


RIESE) | (C63 


M,—-+— 
ee: a ax 


0 as Or 
ax ae a2 =0 (16) 


The result of the subtraction between equations (15) and (16) is again the wave 


equation: 


with solution: 


r(x,2)=f{x+z)+g(x-z) 


OF: 


= fix(i+—_)+a2] + aloes 
Hees) = fills ead) + gle-5p)-af 17) 


B. DEVELOPMENT OF STATIC PRESSURE WAVEFORM 


Assuming an isentropic process, the relation between local pressure and local density 
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Using the previous expression for the local density: 


P=Pol1+s(x,)] 


the local pressure can be written: 


P=P[l+s(x,)]) - PP.[1+ks(x,1)] 


Or: 


P=P,[1+(x,2)] 


where: 


= foils + oe pels ae 
Px!) = FIx(1 M. aot} + G[x(1 he a,f} 


which can be modified to include an harmonic form as: 


~ RoE xy, 1 yO). Gi @hzq_-1) 
p(x,t) areas m, = ara Mm 7} 8) 


Next, we shall assume that the rotating shutter vanes superimpose a harmonic 


pressure perturbation upon the free stream of the form: 


Sy, 


P=P,[1+h coswt] 


where h<1 , at x=0, So: 


p(0,j= F(wt) + G(wt) = h cosut 


The static pressure harmonic perturbations at an upstream distance (x=-L) from the 


source remain consistent with the form of equation (18) and become: 


als Fy} +-cost 24 hi (eee ay 
M, 


P(-L,H= 7 eos 2 =e (1 = M, L 











ae 
P(-Lat)=hcos( S™oost aap, Mest cos( Sint sine 
or 
P(-Lat)=h cos( = =I c0s( op easat + sin( yn | 
or: 


p(-L,t)=h cos(24) cos(—24 
a) aM, 


-wft) 





0 
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with the final expression for the static pressure at a distance L upstream due to the 


pressure oscillation introduced by the rotating shutter vanes: 


P= P[1 +h cos( 2%) cose 4 -wn)] 
eo 0 


OT: 


P=P,[1 +h cos2 cos(- 2-2) 
0 


and the general expression for the static pressure at any distance x and at a given time t, 


due to a harmonic perturbation source at x=0, becomes: 


P(x,t) = P,[ 1 + h cos—cos(—=-wt] (19) 
= Up 


C. DEVELOPMENT OF VELOCITY WAVEFORM 


The above general expression for the static pressure yields: 


p(x, =h cos——~cos(—— - w!) 
Uy 


and since: 


P=P,[1 + ks(z,0) 


s(x,t) has the following form: 
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s(x,t) = ‘ cos——cos( == - cat) (20) 


Differentiation of equation (20) with respect to t, yields: 


Os(x,t) _— ho OX OX_ 71 
= : aera Ts wt] (21) 


while differentiation, again of equation (20) now with respect to x, yields: 


EN UO) a OPS ees ay) Oe ae Ee ra (22) 
ax ka, Uo KU, % U, 


Substitution of equations (21) and (22) into equation (9), yields: 


Or(x,f) _ ho Ba i 


ax kay a Uy 





Integrating the above expression of r(x,t), with respect to x: 


r(x,t) = = f{sin Scos( F*-wn)ds + CW 
or: 
r(x.) = o i sin 2 FE on) osin( SF +o) + CO) 
or: 
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hw fy. | 
ret)= if eink F- (Mla 1)= ct] oxink (Ag ~1) + ori + Ci) 
and finally: 


7 M 
r(x,t) = ou ee 
2k Myt+1 





cos[~(M,+1)-w#] + 
Uy 





+ 


M, 
oA cost My D+we }+ C(t (23) 


Now the rotating shutter vanes superimpose on the free stream a harmonic pressure 


perturbation of the form: 


U = U,{1 + ¢ cos(wt)] 


where q<1 , at x=0 , So: 








r(0,1) = q cos(wt) 
and: 
r(0,t) = at [ mc cos(-wt) + mc cos(wt) J + C(t) 
2k =M,+1 M,-1 
thus: 
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=) eas 





C(t) = (q+ 4 
o 


Substitution of the above expression for C(t) into equation (23), yields the final 


expression for r(x,t): 








ras) = A o*1)-af] + 
Me-1 
+ (M,+ 1)oos[ [= (Mp1) +0 + 2p OO s Mgleasteo } (24) 
0 0 


and the final expression for the velocity in the test section at any distance x from the 


rotating shutter valve and at any time t, becomes: 





h ™M 
u = Uy - Uy = { (M,- TT +1)-of] + 
o7 
gk(Mo-1) 
+ (Mg* Ieost T(Mo-1)+ of + —— (25) 


Now for free stream Mach numbers below 0.1: 


M,<01 - Mo<00l1 - M-1~-1 


equation (25) can be written: 


= U, + U, aye Ue “Neos F(My+1)- wt] + 





+ (M 0* Neosl (Ma -1)+wt] + ls ey ay } (26) 


AM, 
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VI. EXPERIMENTAL PROCEDURES IN THE INVESTIGATION OF THE 


OSCILLATING FLOW FIELD IN THE CLEAR-TUNNEL TEST SECTION 


In order to measure the static pressure oscillations at the same time as the velocity 
oscillation in the clear tunnel test section, a Static pressure probe and a hot wire probe were 
used. Both probes were installed at the same distance of 10.2 feet from the rotating shutter 
vanes and two inches apart, almost in the middle of the test section, just in front of the 
circulation control model future location. 

Exactly the same instrumentation setup designed for the circulation control circular 
cylinder was used to obtain the experimental data. Thus the static pressure probe was 
connected to port number one of the Scanivalve® with the same tubing configuration used 


for its dynamic calibration. 


A. CALIBRATION PROCEDURES 
1. Pressure Measurement System Calibration 
In addition to the test section static pressure the Scanivalve® received 
atmospheric pressure for calibration purposes. With the zero port of the Scanivalve® unit 
Open to the atmosphere, the output of the transducer was adjusted to be equal to zero 
Volts. Then with application of an external static pressure input of 10 inches of water (52 
psf) above atmospheric pressure, the pressure transducer voltage gain was adjusted to read 


1.5 Volts. Thus the voltage output of the Scanivalve® would represent: 


46 


The maximum voltage span of the Scanivalve® was set at plus or minus 5 Volt 
thus pressures equal to P,,, plus or minus 173.33 psf (1.2 psi) were able to be measured 
without over-voltage of the pressure acquisition system. The SRA 1200™ recorder board 
(Channel 1) was also calibrated for a true zero reading at zero Volts input, and set to 
accept plus or minus 2046 digital counts at plus or minus 5 Volts. 

2. Velocity Measurement System Calibration 

The hot-wire anemometer was adjusted at zero velocity, by placing a cover over 
the hot-wire probe, and at maximum gain to obtain an output of Zero Volts. Then with the 
wind tunnel running at 100 ft/sec (the flow was adjusted to generate a water column 
corresponding to 100 ft/sec in the water micromanometer, based on local barometric 
pressure and temperature) and after selecting the 5 Volts scale, the hot-wire output was 
adjusted to obtain a reading of 0.8 full-scale Volts using the gain control knob, 
corresponding to an actual output of 4.0 Volts. 

The SRA 1200™ (Channel 2) recorder board once again was calibrated for a true 
Zero reading at zero volts input and set to accept plus or minus 2046 digital counts at plus 
or minus 5 volts. Thus to avoid over-voltage of the data acquisition system, a maximum 


of 5 Volts, corresponding to 125 feet per second, was used for velocity measurements. 


B. DATA ACQUISITION 

The pressure and velocity signals were output continuously from the Scanivalve® 
pressure transducer and the hot-wire anemometer, respectively, as analog voltages 
corresponding to a percent deflection of full scale, where full scale refered to a full scale 


deflection of the recorder boards. The time-dependent experimental data were recorded 
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and stored in order of occurrence on the 64-K capacity recorder boards. A data acquisition 
program, written in Quick Basic”, was cued by the electric signal which was tripped 
mechanically by the rotating shutter vanes, and armed the two recorder boards. When the 
boards were triggered, 2048 samples of pressure and velocity were recorded 
simultaneously over a period of 1024 milliseconds. This sample rate was programmed to 
match oscillation frequencies in a range between 10 and 50 Hz. 

By calibrating the pressure signal recording board response for a maximum signal 
of plus or minus 173.33 psf (1.2 psi), a capability to resolve 0.085 psf (0.0006 psi) was 
achieved. The velocity recording board was calibrated for a maximum velocity of plus or 
minus 125 ft/sec; therefore its resolution was 0.061 ft/sec. 

The signals were averaged to permit a separation of the time dependent portion of 
each signal from its respective mean value. This resulted in a discrete representation of 
both the oscillating pressure and. the corresponding velocity over a few cycles of flow 
oscillations. A complete program listing appears in Appendix B. 

The operator supplied information to the program regarding tunnel operating 
conditions and ambient conditions for inclusion in the output data file header. The tunnel 
velocity, pressure and velocity signals were monitored continuously to ensure proper 
Operation throughout the run. For steady flow operation, after the rotating shutter vanes 
were secured fully open in the horizontal position, a flow speed of 100 ft/sec was 
established as measured by the Pitot-Static tube and the water micromanometer. Oscillating 
flow fields were investigated at a range of oscillation frequencies of 14 to 40 Hz at 


increments of two Hz. 
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VII. DATA ANALYSIS 


A. DATA REDUCTION 

In this investigation the rotating shutter vanes included three-inch wide blades, 
producing 50 percent closure. After conducting preliminary tests, it became evident that 
the mean value of the free stream velocity and the mean value of the static pressure 
measured from the Scanivalve® did not change with the flow oscillating frequency. The 
objective of these tests was accordingly expanded to include correlation of the 
experimental data with flow parameter variation. The value of 80 ft/sec was measured for 
the mean velocity at a frequency range of 15 to 40 Hz, and the value of 9.9 psf for the 
mean Static pressure was measured by the Scanivalve® for the same frequency range. 

After the time dependent portions of the signals were separated from the mean 
values, unsteady hot-wire velocity signals showed the presence of high frequency noise 
attributable to various sources such as the blower fan. A digital low-pass filter was 
included in the digital processing of the velocity signal with an equivalent low-pass filter 
time constant of five times the sampling frequency, with the algorithm set for a sampling 
rate of about 2000 samples per second. The pressure signal was filtered by the Scanivalve® 
Signal conditioning amplifier and for this reason a clean signal was obtained. 

The pressure data were corrected for phase lag and amplitude attenuation to account 
for the static pressure probe configuration, using the calibration curves for the probe from 
Chapter V. The unsteady pressure and velocity data were compared together versus time, 
using the same time and amplitude scales for each data set, for both phase angle and 


amplitude. 
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B. CORRELATION BETWEEN ANALYTICAL AND EXPERIMENTAL RESULTS 

Velocity and static pressure waveforms showed sinusoidal behavior except around 
25 Hz where a resonance occurred, deforming the velocity wave. It was estimated that the 
resonance appeared due to wave reflection on the constant pressure boundary at the tunnel 
inlet (organ pipe effect). 

After the waveforms for velocity oscillations were plotted, it was apparent that the 
amplitude of oscillation did not changed significantly with the oscillation frequency. With 
some small variations the velocity amplitude maintained the same value of 4 ft/sec 
throughout the whole frequency spectrum that was tested. This value corresponds to 5 
percent of the mean velocity, which was measured to be 80 ft/sec. But when the corrected 
waveforms for static pressure oscillations were plotted, the amplitude of oscillation 
displayed periodic dependence upon the oscillation frequency. The maximum value for the 
Static pressure amplitude of 7.9 psf was observed at the oscillation frequency of 34 Hz 
corresponding to approximately 80 percent of the mean static pressure measured by the 
Scanivalve® (9.9 psf), while the minimum value of 2.0 psf was observed at 18 Hz. 

The above experimental results were used to estimate the amplitude multiplication 
coefficients q and h in the analytical expressions for velocity and pressure harmonic 


perturbations, superimposed on the free stream from the rotating shutter vanes: 


U = U,( 1+ q cosot ) : q<1 


and: 


P=P,.(l+hcosot) , Ah<l 
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Thus the maximum percentage variation of static pressure amplitude was used to 
approximate the amplitude multiplication coefficient h , and the constant percentage 
variation of the velocity amplitude was used to approximate the amplitude multiplication 
coefficient q. Changes to the above coefficients were observed when the flow blockage, 
introduced by the rotating shutter vanes, was increased by increasing the width of the 
blades. 

The analytical expressions developed in Chapter V for the velocity and static 
pressure as functions of the oscillation frequency @ , of the distance from the rotating 
Shutter vanes x , and time , were evaluated using the waveforms obtained from the 
corrected experimental data. Because the hot-wire and the static pressure probes were 
installed at the same distance, L=10.2 ft, from the rotating shutter vanes, it was not possible 
to verify the distance dependance of velocity and static pressure oscillations in the test 
section. The terms "local velocity" and "local static pressure” oscillations were more 
appropriate to be used for the experimentally obtained waveforms, denoting only velocity 
and pressure variations in the vicinity of the future model location, approximately 9 ft from 
the perturbation generator, the rotating shutter valve. 

The analytical results for local velocity and local pressure were plotted versus time 
using the same plot format and at exactly the same frequencies (14 Hz to 40 Hz with 
increments of two Hz) for which the experimental investigation was conducted. The 
analytically obtained waveforms were then superimposed on the existing experimental data 
for direct comparison. Figures 15, 16, 17, 18, 19, and 20 illustrate typical waveforms 
observed for the frequencies 14, 20, 24, 30, 34, and 40 Hz respectively. 

Analytical and experimental data appear to have the same phase angle and the same 


amplitude for all six different frequencies. The local velocity amplitude remained almost 
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constant with increasing oscillation frequency at 4 ft/sec while the local pressure amplitude 
started to decrease at 14 Hz, reaching its minimum value around 20 Hz; it then started to 
increase, reaching its maximum value (8 psf) around 34 Hz; and then it started to decrease 
again. Direct comparison between analytical and experimental results was also made for 
the intermediate frequencies verifying once more the above observations. The relevant 
graphs are listed in Appendix A. 

In order to obtain a more detailed picture of amplitude variations with increasing 
oscillation frequency, the amplitude of local velocity and the amplitude of local static 
pressure were analytically calculated for frequencies between 10 Hz and 100 Hz, and then 
plotted as functions of the oscillating frequency. The behavior of both amplitudes is 
illustrated in Figure 21. 

The significant information appears to lie in the periodic nature of these functions. 
The local static pressure amplitude changes continuously with increasing frequency, 
jumping up between 4 and 80 percent of the value almost every 35 Hz. The local velocity 
amplitude demonstrated a behavior like a sinusoidally amplitude-modulated wave with 
a maximum value 9 percent and a minimum value 1 percent of the mean value. But from 
a macroscopic point of view, local velocity amplitude appears to remain constant at 4 
percent of the mean value. In Figure 21 a comparison was attempted between analytic 
results of local static pressure and velocity amplitudes as a function of frequency. 

Local velocity amplitude changes with frequency are very small compared to local 
static pressure amplitude changes and for this reason local velocity amplitude may be 
considered as a constant for a wide range of oscillation frequencies. Figures 22 and 23 
provide a direct comparison between analytical and experimental results for amplitude 


variations. It is apparent that analytical results correlated well with the experimental 
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OSCILLATIONS AT 14 Hz 
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Figure 15 Velocity And Static Pressure Waveforms At 14 Hz 
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OSCILLATIONS AT 20 Hz 
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Figure 16 Velocity And Static Pressure Waveforms At 20 Hz 
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SCILLATION AT 24 Hz 
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Figure 17 Velocity And Static Pressure Waveforms At 24 Hz 
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OSCILLATIONS AT 30 Hz 
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Figure 18 Velocity And Static Pressure Waveforms At 30 Hz 
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OSCILLATIONS AT 34 hz 
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Figure 19 Velocity And Static Pressure Waveforms At 34 Hz 
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OSCILLATIONS AT 40 Hz 
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Figure 20 Velocity And Static Pressure Waveforms At 40 Hz 
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AMPLITUDES vs, OSCILLATION FREQUENCY 
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21 Static Pressure And Velocity Amplitude Variat 
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Figure 22 Experimental Results For Velocity Amplitude 
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Figure 23 Experimental Results For Static Pressure 
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measurements. From this one might conclude that the pressure amplitude at the source 


(the rotating shutter vanes), h , was itself not strongly dependent on frequency. 
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VIII. CONCLUSIONS 

A Circulation Control model representing a generic, two-dimensional form of a 
tailboom configuration for a no tail rotor (NOTAR™) helicopter was designed and 
fabricated for subsequent testing in the low-speed oscillating flow wind tunnel. A 
frequency response calibration was performed of the unsteady pressure data acquisition 
system to determine the phase lag and amplitude reduction introduced from the 
transmission line connecting the pressure orifices on the model surface to the remote 
pressure transducer inside a Scanivalve®, for a variety of configurations. The calibration 
procedure was also performed on a static pressure probe used to obtain the static pressure 
inside the clear-tunnel test section. An existing real-time data acquisition system was used 
during both the instrumentation and wind tunnel calibrations. Modifications to existing 
software were made as needed. 

Analytical expressions were developed for static pressure and velocity oscillations 
in the tunnel test section, introduced by the harmonic perturbations superimposed upon 
the free stream by the rotating shutter vanes. Comparisons were made between the analytic 
predictions and experimental data obtained during the clear tunnel test section calibrations. 
Data included measurements of local velocity using a constant temperature hot-wire and 
local static pressure using the dynamically calibrated static pressure probe. 

From the direct comparison between experimental and analytical results the 


following conclusions may be drawn: 
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1. Analytical predictions for the oscillating flow field inside the tunnel test section 
strongly agreed with the experimentally obtained waveforms for local velocity and local 
Static pressure. 

2. Changing oscillation frequency produced no significant amplitude variations 
of the local velocity perturbations. The local velocity amplitude remained almost constant 
at 5percent of the mean value as applied frequency was varied. 

3. Oscillation frequency induced significant periodic amplitude variations of local 
Static pressure oscillations with a maximum value 80 percent of the mean value appearing 
first at the oscillation frequency of 18 Hz and followed by peak amplitudes thereafter at 
frequency intervals of 35 Hz. 

The presentation of these results is the first indication of velocity and static pressure 
details in the oscillating flow wind tunnel. Follow on research is suggested to investigate 
the tunnel dynamic behavior at a second tunnel station located at a different axial position. 
Correlation of the static pressure and velocity perturbations between two axial locations 
will provide added credibility to the cos(@x/a,) amplitude modulation factor shown on 
equation (19). Furthermore, it is suggested that confirmation of the dynamic behavior also 
be made at a second value of harmonic blockage. 

The understanding of the flow phenomena in the oscillation frequency low-speed 
wind tunnel will improve the quality of the follow on tests on the circulation control 


circular cylinder model. 
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i1PPENDIX A: Velocity and Static Pressure Waveforms for Oscillation Frequencies 


16, 18, 22, 26, 28, 32, 36, and 38 Hz 


Local velocity and static pressure waveforms are presented for direct comparison 
| with analytical predictions for the above intermediate frequencies, frequencies that were 


not listed in Chapter VI. 
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OSCILLATIONS AT 18 Hz 
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OSCILLATIONS AT 22 Hz 
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OSCILLATIONS AT 26 Hz 
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OSCILLATIONS AT 28 Hz 


—s 
— 
is 
== 
row 
a. 
md 
cy 





awn om DH H D0 OH HD M 


0 








etd Pee 
ogg 





naWneH 0 HH 0 H HW M 


0 


QS PS OS BS GS OS ES CO GS Oo 


gz co Ot AS OI =F OF OP = 


ssd) aunssaud 


TIME {msec} 


70 


OSCILLATIONS AT 92 hz 
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OSCILLATIONS AT 96 Hz 
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OSCILLATIONS AT d8Hz 
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APPENDIX B: Data Acquisition / Data Reduction Program for Clear Tunnel Test 


Section Survey 


‘PROGRAM TKDATA This program contains SUBROUTINES to control the SCANIVALVE, 
‘setup the SRA 1200 RECORDER boards, take dat, download that 
‘data into arrays VRAW,, representing the raw velocity data, 
‘PRAW, representing the raw pressure data, process that data, 
‘and write the output to file and/or plot it on the screen. 
‘These SUBROUTINES and their respective functions are: 


‘HOME - Returns SCANIVALVE to HOME 

‘HEADER - Inputs data for a file header 

‘PORTNUM - Returns SCANIVALVE number (PORTNO) 

‘INIT - Initializes recorder boards 

‘RECORD - Records data 

‘getstartadd - Returns starting address of record 

‘getdat - Returns data samples 

‘PLOT - Plots data on screen 

‘LOAD - Loads two arrays, VRAW and PRAW, (in pseudo- 
‘complex format) 

‘MEAN - Computes the mean of the arrays PRAW and VRAW, 
‘subtracts the means and scales the values 

‘PORT - Advances the SCANIVALVE port 

‘ARRAY - Loads mean values of pressure and velocity 
‘into vectors PMEAN and VMEAN and pressure 
‘and velocity data into arrays P and V 

‘DISK - Outputs pressure, velocity and header data 
‘to a file called DATA.XXX where XXX is the 


“Run Number. 

VSCALE = 125! 
PSCALE = 144! 
OPTION BASE 1 ‘Sets lower bound of subscripts = 1 
DIM physicaladdress(2) ‘Physical adresses of the recorders 
DIM inputscale(2) ‘Full scale input voltage for each 

‘channel. May be 1, 2, 5, or 10 volts 
DIM starts(2) ‘Array containing starting addresses 
DIM VRAW(4096) ‘Array containing raw velocity data 
DIM PRAW(4096) ‘Array containing raw pressure data 
DIM P(100) ‘Alray containing pressure data 
DIM V(100) ‘Array containing velocity data 
GOSUB HEADER ‘Calls for header information 
GOSUB HOME ‘Sets SCANIVALVE to PORT 0 
GOSUB PORT ‘Advances SCANIVALVE to PORT 1 
GOSUB INIT ‘Use once before using recorders 


‘RECORDER CONTROL PARAMETERS 
‘These parameters are: 

rate, the sample rate in mico-seconds, from 1 to 32767 
‘samples, the number of samples per record, from 1 to memmax 
‘  psamples, the number of samples that preceed the trigger from 
‘ 1 to memmax, but less than samples. 
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‘  mode§, trigger slope control, either "PLUS" or "MINUS" 

‘ setpoint, for trigger threashold, range is -100 to 100 % 
optional: 

‘  inputscale(n) controls programmable gain for channel n 

: Valid with values of 1, 2, 5 and 10; n may be either 1 or 2. 


SAMPLES = 2048 ‘Number of samples per record, (must 
‘be of the form 2**n for processing) 
psamples = 1000 ‘Number of pretrigger samples 
setpoint = 10 ‘Trigger setpoint in % of full scale 
mode$ = "PLUS" ‘Trigger mode$, (either "PLUS" 
‘or "MINUS") 
rate = 1000 ‘Sample rate in u-seconds per point 
inputscale(1) = 5 ‘full scale for channel 1 is 5 volts 
inputscale(2) = 5 ‘full scale for channel 2 is 5 volts 
FOR J = 1 TO 20000 ‘Delays data acquisition to allow 
NEXT J ‘SCANIVALVE to settle 
GOSUB PORTNUM ‘Reads current port number 
GOSUB RECORD ‘Records data for one port 
GOSUB PLOT ‘Plots recorded data for debugging 
GOSUB LOAD ‘Loads data into arrays VRAW & PRAW 
GOSUB PORT ‘Advances SCANIVALVE to next port 
GOSUB MEAN ‘Computes the mean of each array, 


‘subtracts the means (PBAR and VBAR), 
‘and scales all values with PSCALE 

‘and VSCALE so that output is in psf 
‘and ft/sec. Arrays remain named PRAW 


‘and VRAW 

GOSUB ARRAY ‘Loads PBAR and VBAR into vectors 
‘and PRAW and VRAW into arrays 

GOSUB DISK: “Writes all header data, the vectors 


‘PMEAN and VMEAN and the arrays P and 
‘V to disk ina file RUNXXX.DAT 
END ‘End of application code... 


‘ SUBROUTINES FOR PROGRAM TKDATA 


; 


INIT: ‘Use this routine once at startup for variable and hardware initialization 


k = 1024 ‘kis 2 ™* 10 

SAMPLES = 2000 ‘init # of samples per record 

psamples = INT(SAMPLES / 4) ‘init presamples at .25 total record 
setpoint = 50 ‘trigger setpoint in % of full scale 

mode$ = "PLUS" ‘trigger mode$ is either "PLUS" or "MINUS" or - 
rate = 4 ‘default sample rate in us. 

memmax = (64 * k)- 1 ‘recorder memory size 

tbe = &H205 ‘time base control register address 

cwr = &H21E ‘8253 control write register 

ctrO = &H206 ‘8253 counter #1 

ctrl = &H20E ‘8263 counter #2 

ctr2 = &H216 8253 counter #3 

clkie = 1 ‘for hardware hook not being used 
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FOR J = 1TO2 


‘set up list for physical address starting at &h200 
physicaladdress(J) = &H200 + J * 32) - 32 


‘Read board jumpers and set inputscale(j) 
‘Absent boards get inputscale = 0 
‘Others will return 1, 2,5 or 10 


x = INP(physicaladdress(J) + 5) ‘read jumpers w4, w5, w6, w7 
x = x AND &H70 


SELECT CASE x 


CASE 96 
inputscale(J) = 1 


CASE 80 
inputscale(]) = 2 


CASE 64 
inputscale(]) = 5! 


CASE 48 
inputscale(J) = 10! 


CASE ELSE 
inputscale(J) = 0 


END SELECT 
NEXT J 


‘Initialize timebase and run through one cycle: 
rate = 1000 


‘Toggle reset line 
OUT tbe, & HEF 
OUT tbe, & HE3 


‘Set counter 0 to mode 3 
OUT ewr, &H36 


‘Set sample rate 
A = rate * 2 

ha = INT(A / 256) 
la = A - (ha * 256) 
OUT ctr, la 

OUT ctr, ha 


‘Set counter 1 to mode 3 
‘And divide by 4 


OUT cwr, &H76 


OUT ctr1, 4 
OUT ctr1, 0 
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‘Set counter 2 to mode 5 
OUT cwr, &HBA 


‘Set number of samples 

B = (SAMPLES - psamples) / 4 
hb = INT(B / 256) 

lb = INT(B - (hb * 256)) 


OUT ctr2, lb 
OUT ctr2, hb 


‘Set trigger level 
A = INT(setpoint * 1.28 * -1) + 127-7 
OUT &H203, A 


‘Arm timebase in chosen mode 


‘First create timebase control word 
tbctw = & HF7 


‘clkie always equals 1 
IF clkie = 1 THEN tbctw = tbctw AND &HFB 


‘Test mode$( PLUS or MINUS ) 

IF mode$ = “MINUS” THEN tbctw = tbctw AND &H7F 
OUT tbc, tbctw AND &HED 

‘Wait loop to get pretrigger data 

ptdelay = psamples * rate * .000001 

A = TIMER + ptdelay + 1 


‘This loop generates and error when timing through midnight 
WHILE (TIMER < A AND INKEY$ = “"): WEND 


OUT tbe, thctw AND &HFE 
OUT tbc, tbctw 


‘Force a trigger 

OUT tbe, tbctw OR 8 
FOR J = 1 TO 500: NEXT J 
OUT tbc, tbctw 


‘This ends the time base first cycle and initialization.. 


RETURN 


RECORD. ‘recording sequence control 
CLS 
“VIEW 
"WINDOW (0, 2048)-(511, -2048) 
“LINE (50, 2000)-(461, -1850), lincol, B 


LOCATE 3, 27 


PRINT "RECORDING SEQUENCE CONTROL " 


LOCATE 8, 25 
PRINT "1. ) Initialize Recorder Hardware " 


‘First set the gain on boards 1 and 2. 
FOR J = 1 TO 2 
SELECT CASE inputscale(J) 


CASE 1 
mask = &H9F 


CASE 2 
mask = & HDF 


CASE 5 
mask = &HBF 


CASE 10 
mask = &HFF 


CASE ELSE 
mask = &HFF 


END SELECT 


IF J = 1 THEN mask] = mask 
OUT physicaladdress(J) + 5, mask 


NEXT J 


‘Toggle reset line 
OUT tbe, & HEF AND mask1 
OUT tbe, & HE3 AND mask]! 


‘Set counter 0 to mode 3 
OUT cwr, &H36 


‘Set sample rate 
A = rate * 2 


ha = INT(A / 256) 
la = A - (ha * 256) 
OUT ctr, la 
OUT ctr0, ha 


‘Set counter 1 to mode 3 
‘And divide by 4 


OUT cwr, &H76 
OUT ctrl, 4 
OUT ctrl, 0 


‘Set counter 2 to mode 5 
OUT cwr, &HBA 
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‘Set number of samples 

B = (SAMPLES - psamples) / 4 
hb = INT(B / 256) 

lb = INT(B - (hb * 256)) 


OUT ctr2, Ib 
OUT ctr2, hb 


‘Set trigger level 
A = INT(setpoint * 1.28 * -1) + 127 - 7 
OUT &H203, A 


‘Arm timebase in chosen mode 


‘First create timebase control word 
tbctw = &HF7 AND maskl 


‘clkie always equals 1 
IF clkie = 1 THEN tbctw = tbctw AND &HFB 


“Test mode$( PLUS or MINUS ) 
IF mode$ = "MINUS” THEN tbctw = tbctw AND &H7F 


LOCATE 10, 25 
PRINT "2. ) Pretrigger Delay ” 


OUT tbc, thctw AND &HED 
REM out tbc,tbctw 


‘Wait loop to get pretrigger data 
ptdelay = psamples * rate * .000001 
A = TIMER + ptdelay + 1 


‘This loop generates and error when timing through midnight 
WHILE (TIMER < A AND INKEY$ = ""): WEND 


LOCATE 12, 25 
PRINT "3. ) Arm Recorder...wait for trigger” 


OUT tbe, tbctw AND &HFE 
OUT the, tbctw 


BS = 
A=] 


WHILE (A <> 0 AND INKEY$ = "") 
A = INP(&H205) AND 3 


IF INKEY$ <> "" THEN 
OUT tbc, tbctw OR 8 
FOR J = 1 TO 500: NEXT J 
OUT tbc, tbctw 
END IF 

WEND 


LOCATE 14, 25 
PRINT “4. ) Triggered ...wait until recording complete 
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WHILE (A = 1 AND INKEY$S = ™") 
A = INP(&H205) AND 1 
WEND 


LOCATE 16, 25 
PRINT "6. ) Recording port number..."; PORTNO 
PRINT "6. ) Hit any Key to Exit" 


A$ =" 
WHILE A$ =" 
A$ = INKEY$ 
WEND 


CES 


RETURN 


getstartadd: ‘read 18 bit stop address from 3 ports and convert to number 
‘only 16 bit stop address from 2 ports for 64 k sample memory 


stopadd = 0 


‘bits 16 and 17 

REM sadd1 = INP(physicaladdress(i) + 2) AND &H3 
REM sadd1 = sadd1 * 64 * 1024 

REM stopadd = sadd1 


“bits 8 thru 15 

sadd1 = INP(physicaladdress(I) + 1) 
sadd1 = sadd1 * 256 

stopadd = stopadd + sadd1 


‘bits 0 thru 7 
stopadd = INP(physicaladdress(I)) + stopadd 


‘Calculate starting address by subtracting samples 
starts(I) = stopadd - SAMPLES - 4 


‘Starting address is never negative, so fix it if needed. 
IF starts(I) < 0 THEN starts(I) = starts(T) + 65536 


RETURN 
‘GETDAT and gd1 both get a single value from the recorder memory, 
‘this value is returned in the variable DAT. DAT is a 
‘number that ranges from -2048 to +2047. 
‘GETDAT will initialize the memory pointer with the value in BRDST. 
‘However, this memory pointer need only be initialized one time per 


‘record transfer. Then the calling routine can continously call gd1 
‘and the recorder memory pointer will auto increment after each call 


getdat: 
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IF brdst > memmax THEN brdst = brdst - memmax 
IF brdst < 0 THEN brdst = brdst + memmax 


REM str = INT(brdst / 65536) 

REM brdst = brdst - (str * 64 * 1024) 
strhi = INT(brdst / 256) 

strlo = brdst - (strhi * 256) 


OUT physicaladdress(I), strlo 
OUT physicaladdress(I) + 1, strhi 
REM OUT physicaladdress(i) + 2, str 


gd: 
A = INP(physicaladdress(I) + 3) 
B = INP(physicaladdress(I) + 4) 
B = (B AND &HF) * 256 
DAT = (A + B - 2048) * -1 


RETURN 


PLOT: ‘This routine will plot the arrays to the screen. 


SCREEN 2 ‘2 for cga, 3 for hercules and 
‘9 for ega. (Only required for 


‘plotting) 
WINDOW (0, 2048)-(SAMPLES, -2048) 


‘plot data on screen 
FOR I= 1 TO 2 


GOSUB getstartadd ‘Determines the starting point of 
‘the record in memory 


brdst = starts(I) ‘Initializes memory pointer and gets 
GOSUB getdat ‘first data point 
PSET (0, DAT) 


FOR J = 1 TO SAMPLES 
GOSUB gd1 ‘gets sample in dat and auto inc. ptr 


LINE -J, DAT) 
NEXT J 
NEXT I 


RETURN 


‘ 


LOAD: ‘Load data into an array 
FOR I = 1 TO 2 


GOSUB getstartadd ‘Initializes boards for data 


brdst = starts(I) ‘downloading and gets first 
GOSUB getdat ‘sample 
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FOR J = 1 TO (2 * SAMPLES) STEP 2 


GOSUB gdl ‘Gets sample and increments pointer 
IF I = 1 THEN 
VRAW() = DAT 
VRAW( + 1) = 0! 
END IF 
IF I = 2 THEN 
PRAW() = DAT 
PRAW] + 1) = 0! 
END IF 


NEXT J 
NEXT I 


RETURN 
HEADER: 


‘This subroutine collects run information for a file header and 
‘data processing 


CLS ‘Clear screen 

LOCATE 2, 10 ‘Locates text 

PRINT "RUN DATA INPUT" 

LOCATE 4, 20 

PRINT "1. Input a Run Number consisting of the Model” 
LOCATE 5, 25 


PRINT "Number and the Angle of Attack in the form XYY.” 
INPUT RUNNOS 

LOCATE 6, 25 

PRINT "RUN NUMBER - *; RUNNO$ 

LOCATE 8, 20 

PRINT "2. Input Date” 

INPUT DAIT$ 

LOCATE 9, 25 

PRINT "DATE-"; DAIT$ 

LOCATE 11, 20 

PRINT "3. Input Time” 

INPUT TYME$ 

LOCATE 12, 25 

PRINT "TIME-’; TYME$ 

LOCATE 14, 20 

PRINT "4. Input Model Description” 

INPUT MODELS 

LOCATE 15, 25 

PRINT "MODEL-"; MODELS 

LOCATE 17, 20 

PRINT "5. Input Angle of Attack (+/-) in Degrees” 

INPUT AOA 

LOCATE 18, 25 

PRINT “ANGLE OF ATTACK-"; AOA; “DEGREES” 
LOCATE 20, 20 

PRINT "6. Input Mean Dynamic Pressure in CM of Water” 
INPUT PDYN 

LOCATE 21, 25 

PRINT "MEAN DYNAMIC PRESSURE-"; PDYN; "CM-WATER"” 
LOCATE 23, 20 
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a A 





PRINT "7. Input Indicated Oscillation Frequency (Hz)" 
INPUT OMEGA 

LOCATE 24, 25 

CLS 

LOCATE 3,20 

PRINT "OSCILLATION FREQUENCY-"; OMEGA; "Hz" 
LOCATE 5, 20 

PRINT "8. Input Barometric Pressure (In Hg)" 

INPUT PREF 

LOCATE, 23 

PRINT "BAROMETRIC (REF) PRESSURE-"; PREF; "IN HG" 
LOCATE 9,20 

PRINT "9. INPUT Freestream Temperature (F)" 

INPUT TEMP 

LOCATE 10, 25 

PRINT "FREESTREAM TEMPERATURE-"; TEMP; "F" 


RETURN 


HOME: 
‘This subroutine advances the SCANIVALVE to HOME 


FOR I = 1 TO 50000 


OUT &H2A8, 2 
NEXT I 
OUT &H2A8, 0 
RETURN 
PORTNUM: 


‘This file contains the text of the subprogram PORTNUM in TKDATA 


‘This subprogram reads the BCD output of the SCANIVALVE, 
‘(PORTNO) 


PORTNO = 255 - INP(&H2A9) 
RETURN 


PORT: 


‘This file contains the text of the subprogram PORT in TKDATA 
‘This subprogram advances the SCANIVALVE one port 


FOR I = 1 TO 1000 

OUT &H2A8, 1 

NEXT I 

OUT &H2A8, 0 
RETURN 


MEAN: ‘Determines the mean values of the two arrays PRAW and VRAW, PBAR 
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‘and VBAR, and replaces the arrays with the mean values subtracted 


VBAR = 0! 
PBAR = 0! 


FOR l= 1 TO (2 * SAMPLES) STEP 2 
VBAR = CDBL(VRAW()) * (VSCALE / 2048) + VBAR) 
PBAR = CDBL(PRAW()) * (PSCALE / 2048) + PBAR) 
NEXT I 


VBAR = VBAR / SAMPLES 
PBAR = PBAR / SAMPLES 


FOR I = 1 TO (2 * SAMPLES) STEP 2 
VRAW(]) = VRAW() * (VSCALE / 2048) - VBAR 
PRAW(I) = PRAW(]) * (PSCALE / 2048) - PBAR 
NEXT I 


RETURN 


ARRAY: ‘This subroutine loads PBAR and VBAR into vectors PMEAN and VMEAN, 
‘and PRAW and VRAW into arrays P and V. (It truncates the 1024 
‘records in PRAW and VRAW into 1000 records, performs an ensemble 
‘average over the 10 ensembles and loads a 100 sample array.) 


PMEAN(PORTNO) = PBAR 
VMEAN(PORTNO) = VBAR 


FOR k = 1 TO 100 
PEA = 0! 
VEA = 0! 
FORJ =0TO9 
PEA = PRAW((2 * k)- 1+ * 100)) 
VEA = VRAW((2 * k) - 1 + J * 100)) 
NEXT J 


P(PORTNO, k) = PEA 
V(PORTNO, k) = VEA 
NEXT k 


RETURN 


DISK: ‘This subroutine creates a file A:RUNXXX.DAT, (where XXX is the 
“RUN NUMBER), and writes all HEADER DATA, PMEAN and VMEAM 
‘and the pressure and velocity arrays, P and V to that file. 


CLS 
LOCATE 10, 15 
PRINT “INSERT A FORMATTED DISK IN DRIVE A:" 


LOCATE 15, 15 
PRINT "PRESS ANY KEY TO WRITE RUN DATA TO FILE" 


WHILE A$ =”" 
A$ = INKEYS 
WEND 


CLS 


LOCATE 10, 10 
PRINT "Run data is being written to file ASRUN” + RUNNOS + ".DAT" 


‘Open the file: 


OPEN "A:RUN" + RUNNOS + ".DAT”" FOR OUTPUT AS #2 


‘Write header data: 


WRITE #2, RUNNOS, DAITES, TYMES, MODELS: ‘Strings! 
WRITE #2, AOA: ‘degrees 

WRITE #2, PDYN * 2.048: ‘psf 

WRITE #2, OMEGA: ‘hertz 

WRITE#2, PREF * 70.73: ‘pst 

WRITE #2, TEMP + 459.7: ‘deg Rankine 


“Next write pressure and velocity data: 
WRITE #2, PMEAN 
WRITE #2, VMEAN 


CLS 


BEEP 


FOR J = 1 TO 100 

WRITE #2, PQ): ‘Pressure array 
NEXT J 
FOR J = 1 TO 100 

WRITE #2, VQ): ‘Velocity array 
NEXT J 


LOCATE 10, 10 
PRINT" = RUN COMPLETE AND DATA WRITTEN TO FILE” 


RETURN 
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APPENDIX C: Data Acquisition / Data Reduction Program for Unsteady Pressure Measurement System 
Calibrations 


‘PROGRAM DYNARES.BAS (DYNAmic RESponse) 


‘The program records the pressure signals from the scanivalve transducer \ 
‘( Channel 2 ) and the reference pressure transducer ( Channel 1 ). 

‘The pressure signals are output continuously from the transducers as DC 
‘voltages corresponding to a percent deflection of full scale of the recorder \ 
‘poard. In both channels, full scale deflection is limited to +/- 5 Volts. | 
‘The program arms the boards and when is triggered ( manually from the 
‘rotating shutter vanes switch ), 1024 samples of pressure are recorded 
‘simultaneously from each channel over a period of 1.024 seconds ( sample 
‘rate 1000 samples/sec ) which spans 10 cycles of pressure oscillation. 

‘The program uses the following subroutines: 


‘HEADER - Inputs data for a file header 

‘INIT - Initializes recorder boards 

“RECORD - Records data 

‘getstartadd - Returns starting address of record 
‘getdat - Returns data samples 

‘PLOT - Plots data on screen 

‘LOAD - Loads two arrays, PIRAW and P2RAW 


‘DISK - Outputs PIRAW, P2RAW and header data 
: to a file called SCAN.XXX where XXX is the scanivalve port number 





P1ISCALE = 144! 

P2SCALE = 144! 

OPTION BASE 1 ‘Sets lower bound of subscripts = 1 

DIM physicaladdress(2) ‘Physical adresses of the recorders 1 

DIM inputscale(2) ‘Full scale input voltage for each 1 

‘channel. May be 1, 2, 5, or 10 volts 

DIM starts(2) ‘Array containing starting addresses 

DIM P1RAW(2048) ‘Array containing raw pressure1 data 

DIM P2RAW(2048) ‘Array containing raw pressure2 data Mi, 
iy 

GOSUB HEADER ‘Calls for header information 7 

GOSUB INIT ‘Use once before using recorders 


f 


‘RECORDER BOARDS CONTROL PARAMETERS 


SAMPLES = 1024 ‘Number of samples per record, (must 
‘be of the form 2**n for processing) 
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psamples = 1000 ‘Number of pretrigger samples 


setpoint = 10 ‘Trigger setpoint in % of full scale 
mode$ = "PLUS" ‘Trigger mode$, (either "PLUS" 
‘or "MINUS") 
rate = 1000 ‘Sample rate in u-seconds per point 
inputscale(1) = 10 ‘full scale for channel 1 is 10 volts 
inputscale(2) = 10 ‘full scale for channel 2 is 10 volts 
FOR J = 1 TO 20000 ‘Delays data acquisition to allow 
NEXT J ‘scanivalve to settle 
GOSUB RECORD ‘Records data for one port 
GOSUB PLOT ‘Plots recorded data for debugging 
GOSUB LOAD ‘Loads data into arrays PIRAW & P2RAW 
‘of length 2048 
GOSUB PHASE ‘Calculates the phase lag between 


‘arrays P2ARAW and P2RAW 


‘GOSUB DISK: ‘Writes all header data, PIRAW and 
‘P2RAW to disk ina file SCNXXX.DAT 


END ‘End of application code...... 


ee ce cc ce ce cc ce ce ce ce ce ce ge ce ee ee ce ce ee ce we ee ce ee ee ee ee ee ee ee ee oe oe a 


‘ SUBROUTINES WHICH CALLED BY PROGRAM DYNARES 


INIT: ‘Use this routine once at startup for variable and hardware init. 


k = 1024 ‘kis 2 ** 10 

SAMPLES = 2000 ‘init # of samples per record 

psamples = INT(SAMPLES / 4) ‘init presamples at .25 total record 
setpoint = 50 ‘trigger setpoint in % of full scale 

mode$ = "PLUS” ‘trigger mode$ is either "PLUS" or "MINUS" or - 
rate = 4 ‘default sample rate in us. 

memmax = (64 * k) - 1 ‘recorder memory size 

tbe = &H205 ‘time base contro] register address 

cwr = &H21E ‘8253 control write register 

ctrO = &H206 8253 counter #1 

ctrl! = &H20E ‘8263 counter #2 

ctr2 = &H216 8253 counter #3 

clkie = 1 ‘for hardware hook not being used 

FOR J = 1 TO2 


‘set up list for physical address starting at &h200 
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physicaladdress(J]) = &H200 + (J * 32) - 32 


‘Read board jumpers and set inputscale(j) 
‘Absent boards get inputscale = 0 
‘Others will return 1, 2, 5 or 10 


x = INP(physicaladdress(J) + 5) ‘read jumpers w4, w5, w6, w7 
x = x AND &H70 


SELECT CASE x 


CASE 96 
inputscale(J) = 1 


CASE 80 
inputscaleJ) = 2 


CASE 64 
inputscale(Q) = 5! 


CASE 48 
inputscale(J) = 10! 


CASE EESE 
inputscaleJ) = 0 


END SELECT 


NEXT J 


‘Initialize timebase and min through one cycle: 
rate = 1000 
‘Toggle reset line 


OUT tbe, & HEF 
OUT tbe, &HE3 


‘Set counter 0 to mode 3 
OUT cwr, &H36 


‘Set sample rate 

A = rate * 2 

ha = INT(A / 256) 
la = A - (ha * 256) 
OUT ctr0, la 

OUT ctr0, ha 


‘Set counter 1 to mode 3 
‘And divide by 4 


OUT cwr, &H76 
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OUT ctrl, 4 
OUT ctr1, 0 


‘Set counter 2 to mode 5 
OUT cwr, & HBA 


‘Set number of samples 

B = (SAMPLES - psamples) / 4 
hb = INT(B / 256) 

lb = INT(B - (hb * 256)) 


OUT ctr2, lb 
OUT ctr2, hb 


‘Set trigger level 
A = INT(setpoint * 1.28 * -1) + 127 - 7 
OUT &H203, A 


‘Arm timebase in chosen mode 


‘First create timebase control word 
tbctw = &HF7 


‘clkie always equals 1 
IF clkie = 1 THEN tbctw = tbctw AND &HFB 


‘Test mode$( PLUS or MINUS ) 

IF mode$ = "MINUS" THEN tbctw = tbctw AND &H7F 
OUT tbc, tbctw AND &HED 

‘Wait loop to get pretrigger data 

ptdelay = psamples * rate * .000001 

A = TIMER + ptdelay + 1 


‘This loop generates and error when timing through midnight 
WHILE (TIMER < A AND INKEY$ = ""): WEND 


OUT tbc, tbctw AND &HFE 
OUT tbc, tbctw 


‘Force a trigger 

OUT tbc, tbctw OR 8 

FOR J = 1 TO 500: NEXT J 

OUT tbc, tbctw 

‘This ends the time base first cycle and initialization. 
RETURN 
RECORD. ‘recording sequence control 

CLS 


“VIEW 
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"WINDOW (0, 2048)-(511, -2048) 
‘LINE (50, 2000)-(461, -1850), lincol, B 


LOGATE 3727 
PRINT “RECORDING SEQUENCE CONTROL " 


LOCATE 8, 25 
PRINT "1. ) Initialize Recorder Hardware " 


‘First set the gain on boards 1 and 2. 
FOR J = 1 TO 2 
SELECT CASE inputscale() 


CASE 1 
mask = &H9F 


CASE 2 
mask = &HDF 


CASE 5 
mask = &HBF 


CASE 10 
mask = &HFF 


CASE ELSE 
mask = &HFF 


END SELECT 


IF J = 1 THEN mask! = mask 
OUT physicaladdress(J) + 5, mask 


NEXT J 


‘Toggle reset line 
OUT tbe, & HEF AND mask1 
OUT tbe, & HE3 AND mask! 


‘Set counter 0 to mode 3 
OUT cwr, &H36 


‘Set sample rate 
A=rate*2 


ha = INT(A / 256) 
la = A - (ha * 256) 
OUT ctr, la 
OUT ctr0, ha 


‘Set counter 1 to mode 3 
‘And divide by 4 
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OUT cwr, &H76 
OUT ctr1, 4 
OUT ctr1, 0 


‘Set counter 2 to mode 5 
OUT cwr, &HBA 


‘Set number of samples 

B = (SAMPLES - psamples) / 4 
hb = INT(B / 256) 

lb = INT(B - (hb * 256)) 


OUT ctr2, Ib 
OUT ctr2, hb 


‘Set trigger level 
A = INT(setpoint * 1.28 * -1) + 127 - 7 
OUT &H203, A 


‘Arm timebase in chosen mode 


‘First create timebase control word 
tbctw = &HF7 AND mask] 


‘clkie always equals 1 
IF clkie = 1 THEN tbetw = tbctw AND &HFB 


‘Test modeS( PLUS or MINUS ) 

IF mode$ = "MINUS" THEN tbctw = tbctw AND &H7F 
LOCATE 10, 25 

PRINT "2. ) Pretrigger Delay ” 

OUT tbc, thctw AND &HED 

REM out tbc,tbctw 


‘Wait loop to get pretrigger data 
ptdelay = psamples * rate * .000001 
A = TIMER + ptdelay + 1 


‘This loop generates and error when timing through midnight 
WHILE (TIMER < A AND INKEY$ = ""): WEND 

E@CATE 12,25 

PRINT "3. ) Arm Recorder...wait for trigger” 


OUT tbc, thctw AND &HFE 
OUT tbc, tbctw 


BS = "" 
A=1 
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WHILE (A <> 0 AND INKEY$ = “") 
A = INP(&H205) AND 3 


IF INKEY$ <> “" THEN 
OUT tbe, tbctw OR 8 
FOR J = 1 TO 500: NEXT J 
OUT tbe, tbctw 
END IF 

WEND 


LOCATE 14, 25 
PRINT "4. ) Triggered ...wait until recording complete” 


WHILE (A = 1 AND INKEYS = “") 
A = INP(&H205) AND 1 
WEND 


RETURN 


’ 


getstartadd: ‘read 18 bit stop address from 3 ports and convert to number 
‘only 16 bit stop address from 2 ports for 64 k sample memory 


stopadd = 0 


‘bits 16 and 17 

REM sadd1 = INP(physicaladdress(i) + 2) AND &H3 
REM sadd1 = sadd1 * 64 * 1024 

REM stopadd = sadd1l 


‘bits 8 thru 15 

sadd1 = INP(physicaladdress(I) + 1) 
sadd1 = sadd1 * 256 

stopadd = stopadd + sadd1 


‘bits 0 thru 7 
stopadd = INP(physicaladdress(I)) + stopadd 


‘Calculate starting address by subtracting samples 
Starts(I) = stopadd - SAMPLES - 4 


‘Starting address is never negative, so fix it if needed. 
IF starts(I) < 0 THEN starts(I) = starts(I) + 65536 


RETURN 


‘GETDAT and gd1 both get a single value from the recorder memory, 
‘this value is returned in the variable DAT. DAT is a 
‘number that ranges from -2048 to +2047. 


‘GETDAT will initialize the memory pointer with the value in BRDST. 
‘However, this memory pointer need only be initialized one time per 
‘record transfer. Then the calling routine can continously call gd1 
‘and the recorder memory pointer will auto increment after each call. 


92 














getdat: 
IF brdst > memmax THEN brdst = brdst - memmax 
IF brdst < 0 THEN brdst = brdst + memmax 


REM str = INT(brdst / 65536) 

REM brdst = brdst - (str * 64 * 1024) 
strhi = INT(brdst / 256) 

strlo = brdst - (strhi * 256) 


OUT physicaladdress(I), strlo 
OUT physicaladdress(I) + 1, strhi 
REM OUT physicaladdress(i) + 2, str 


gd: 
A = INP(physicaladdress(I) + 3) 
B = INP(physicaladdress(I) + 4) 
B = (B AND &HF) * 256 
DAT = (A + B - 2048) * -1 


RETURN 


s 


PLOT: ‘This routine will plot the arrays to the screen. 


SCREEN 2 ‘2 for cga, 3 for hercules and 
‘9 for ega. (Only required for 


‘plotting) 
WINDOW (0, 2048)-(SAMPLES, -2048) 


‘plot data on screen 
FOR I = 1TO2 


GOSUB getstartadd ‘Determines the starting point of 
‘the record in memory 
brdst = starts(1) ‘Initializes memory pointer and gets 
GOSUB getdat ‘first data point 
PSET (0, DAT) 


FOR J = 1 TO SAMPLES 
GOSUB gdl ‘gets sample in dat and auto inc. ptr 


LINE -J, DAT) 


NEXT J 
NEXT I 
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LOAD: ‘Load data into an array 


FOR 1=1TO2 
GOSUB getstartadd ‘Initializes boards for data 
brdst = starts(I) ‘downloading and gets first 
GOSUB getdat ‘sample 


FOR J = 1 TO (2 * SAMPLES) STEP 2 


GOSUB gd ‘Gets sample and increments pointer 
IF l= 1 THEN 
PIRAW() = DAT 
PIRAW( + 1) = 0! 
END IF 
IF ] = 2 THEN 
P2ZRAWQ) = DAT 
P2RAWG + 1) = 0! 
END IF 


HEADER: 
‘This subroutine collects run information 
CLS ‘Clear screen 


LOCATE 2, 10 ‘Locates text 
PRINT "RUN DATA INPUT” 


LOCATE 4, 20 

PRINT “Number of run and the port in the form XYY." 
INPUT RUNNO$ 

LOCATE 6, 25 

PRINT "RUN NUMBER - "; RUNNO$ 


LOCATE 8, 20 

PRINT "7. Input Indicated Oscillation Frequency (Hz)" 
INPUT OMEGA 

LOCATE 10, 25 

PRINT “OSCILLATION FREQUENCY-"; OMEGA; "Hz" 


RETURN 


DISK: ‘This subroutine creates a file A:SCNXXX.DAT, (where XXX is the 
“RUN NUMBER AND SCANIVALVE PORT), and writes all HEADER DATA, 
‘PIRAW and P2RAW to that file. 


CLS 


LOCATE 10, 15 
PRINT “INSERT A FORMATTED DISK IN DRIVE A:" 
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LOCATE 15,15 
PRINT "PRESS ANY KEY TO WRITE RUN DATA TO FILE" 


AS =“ 
WHILE A$ =" 
AS = INKEYS 
WEND 


CLS 
LOCATE 10, 10 
PRINT "Run data is being written to file ASSCN” + RUNNO$ + “.DAT" 


‘Open the file: 
OPEN "A:SCN" + RUNNOS + ".DAT" FOR OUTPUT AS #2 
‘Write header data: 


WRITE #2, RUNNOS, DAITES: ‘Strings! 
WRITE #2, OMEGA: ‘hertz 


‘Next write pressure data: 


FOR I = 1 TO (2 * SAMPLES) STEP 2 
WRITE #2, PIRAW()): 
NEXT I 


FOR I = 1 TO (2 * SAMPLES) STEP 2 
WRITE #2, P2ARAW(I): 
NEXT I 


es 


BEEP 
LOCATE 10, 10 
PRINT" RUN COMPLETE AND DATA WRITTEN TO FILE" 


RETURN 


‘ 


PHASE: ‘This subroutine calculates the phase lag in degrees between 
‘the pressure signals from Channel Nol and Channel No2. 


SUMP1 = 0 

SUMP2 = 0 

SP1P2 = 0 

FOR I = 1 TO (2 * SAMPLES) STEP 2 
SUMP1 = SUMP1 + PIRAW(I) * PIRAW(I) 
SUMP2 = SUMP2 + P2RAW(I) * P2RAW()) 
SP1P2 = SP1P2 + PIRAW(I) * P2RAW(1) 

NEXT I 


N = SAMPLES 


P1 = SQR(2 * SUMP] / N) 
P2 = SQR(2 * SUMP2 / N) 
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COSDFl = (2 * SViir2)7 (N= Pi * P2) 
ARCCOSDFI = 3.141593 / 2 - ATN(COSDFI / SQR(-COSDFI * COSDF + 1)) 
DFI = ARCCOSDFI * 180 / 3.141593 


PRINT "THE PHASE LAG BETWEEN CHANNEL 1 AND CHANNEL 2 IN DEGREES IS :" 
PRINT DFI 


RETURN 


‘ 
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